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By using only a fraction of the model space extended beyond current no-core shell-model limits 
and a schematic effective many-nucleon interaction, we gain additional insight within a symmetry- 
guided shell-model framework, into the many-body dynamics that gives rise to the ground state 
rotational band together with phenomena tied to alpha-clustering substructures in the low-lying 
states in 12 C, and in particular, the challenging Hoyle state and its first 2 + excitation. For these 
states, we offer a novel perspective emerging out of no-core shell-model considerations, including a 
discussion of associated nuclear shapes and matter radii. This, in turn, provides guidance for ab 
initio shell models by informing key features of nuclear structure and the interaction. 



Our present-day knowledge of various phenomena of 
astrophysical significance, such as nucleosynthesis, the 
evolution of primordial stars in the Universe, and X-ray 
bursts depends on reaction rates for the stellar triple- 
a process, which can considerably affect, e.g., results of 
core-collapse supernovae simulations and stellar evolu- 
tion models, predictions regarding X-ray bursts, as well 
as estimates of carbon production in asymptotic giant 
branch (AGB) stars pQ. These rates, in turn, are greatly 
influenced by accurate measurements and theoretical pre- 
dictions of several important low- lying states in 12 C, in- 
cluding the second (Hoyle) state and its 2 + excitation 
that continues to foster debate in experimental studies 
[2H6]. Further challenges relate to the long-recognized 
alpha-cluster substructure of these states that has been 
limited to explorations by cluster-tailored [TJ |H] or self- 
consistent [9] microscopic methods. Indeed, these highly 
deformed patterns have hitherto precluded an accurate 
- from first principles - no core shell-model (NCSM) de- 
scription of the Hoyle state and its rotational band [TO] , 
and, only recently, first ab initio state-of-the-art calcu- 
lations have been attempted using effective field theory 
(EFT) and a coarse lattice [TT] , 

In this letter, we report on a first study of these phe- 
nomena in 12 C within a no-core shell- model framework 
with essentially no limitation on the number of harmonic 
oscillator (HO) shells included in the model space. While 
such model spaces remain inaccessible by ab initio shell 
models, we are able to address a long-standing challenge 
[T2] , namely, understanding highly-deformed spatial con- 
figurations from a shell-model perspective. Furthermore, 
this study allows one to gain further insight into the 
many-body dynamics, including the physically relevant 
shape deformations and particle- hole configurations, that 
gives rise to the ground state (g.st.) rotational band (the 



lowest + , 2 + and 4 + states) together with low- lying 
states suggested to have a cluster structure (Oj Hoyle 
state and its 2 + excitation), as well as a third low- lying 
O3" state in 12 C. This analysis greatly benefits from the 
use of a schematic, but fully microscopic, effective many- 
nucleon interaction and is based on the close agreement of 
model outcomes with experiment - including excitation 
energies and other observables such as mass rms radii, 
electric quadrupole moments and E2 transition rates - 
as well as to wavefunctions obtained by ab initio no-core 
shell model (NCSM) calculations [13^ with the JISP16 
realistic nucleon-nucleon (AW) interaction [T4] , 

Above all, the qualitative results of the present study 
provide guidance for ab initio shell model approaches by 
informing key features of nuclear structure and the inter- 
action. In particular, for a description of alpha-cluster 
phenomena, the outcome points to the need for simple 
effective interactions beyond two-body realistic ones, and 
manifests a clear indication that achieving ab initio de- 
scriptions is within the reach of the NCSM. The latter, 
in turn, will bring forward an accurate reproduction and 
reliable prediction of energy spectra and associated tran- 
sition rates that majorly impact astrophysical studies. 

Symmetry-adapted shell-model framework. — 

We employ the no-core symplectic model (NCSpM) for 
symmetry-preserving interactions with Sp(3,R) the un- 
derpinning symmetry [15]. This symmetry is indeed 
found inherent to nuclear dynamics - a result we have 
demonstrated in an analysis of large-scale ab initio 
NCSM applications for 12 C and 16 [16 . The model 
offers a microscopic description of A nucleons in terms 
of mixed shape deformations, directly related to particle 
relative (with respect to the center of mass, CM) position 
and momentum coordinates, r^ and p^, with i = 1, . . . , A. 
It has been successfully applied to 166 Er using Davidson 
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FIG. 2: NCSpM energy spectrum of (a) 12 C, not scaled 
(dotted) and rescaled (solid), and (b) 24 Mg (not scaled). Ex- 
perimental energies for all states are from [29] , except the 
latest results for the 12 C 0^ and 2 J [4] (shaded area shows 
the 2 2 energy range from [2-6 ). B(E2) transition rates are 
in W.u. 



FIG. 1: Sp(3, R) irreps (slices) that comprise the model space 
used for the 12 C NCSpM calculations. Basis states (A/i) of a 
slice are built by 2hQ lp-lh monopole or quadrupole excita- 
tions (set II) over a bandhead. The symplectic bandhead (set 
I) is a SU(3)-coupled many-body state with a given nucleon 
distribution over the HO shells. The corresponding HO en- 
ergy of this nucleon configuration together with the bandhead 
shape, (Ao- /i<x), serve to label the symplectic irrep. 

potential [17J and is a microscopic realization of the Bohr- 
Mottelson collective model [18] , as well as a multiple HO 
shell generalization of Elliott's SU(3) model [T9] . 

The NCSpM utilizes a symplectic basis (for details, see 
[20]), which is related - via a unitary transformation - to 
the three-dimensional HO (m-scheme) many-body basis 
used in the NCSM [21]. The conventional NCSM ba- 
sis spaces are constructed using HO single-particle states 
and are characterized by the Ml oscillator strength as well 
as by the cutoff in total oscillator quanta, iV max , above 
the lowest energy configuration for a given nucleus. In- 
deed, the NCSpM employed within a full model space up 
through iV max , will coincide with the NCSM for the same 
iV max cutoff. It is therefore clear that the present study, 
while down-selecting to the most relevant configurations, 
provides the first shell-model calculations carried beyond 
current NCSM limits, namely, up through 7V max = 20, 
the model space we found sufficient for the convergence 
of results. These important configurations are chosen 
among all possible symplectic Sp(3,R) irreducible repre- 
sentations (irreps) within the model space. 

The Sp(3, R) irreps divide the space into 'vertical slices' 
that are comprised of basis states of a definite shape 
specified by the deformation (A /i) quantum numbers of 
SU(3), which is embedded in Sp(3,R) (Fig. 0. The 
simplest cases, namely, (0 0), (AO), and (0/i), describe a 
spherical, prolate, and oblate shape, respectively, while a 
general nuclear state is typically a superposition of sev- 



eral hundred various triaxial shapes. The basis states 
are built over a bandhead (Fig. [I] Set I) by consecutive 
2MI lp-lh (1-particle-l-hole) excitations (Fig. [I] Set II), 
together with a smaller 2MI 2p-2h correction for elimi- 
nating the spurious CM motion (not shown in the figure). 
In NCSpM, to eliminate the spurious CM motion, we use 
symplectic generators constructed in relative coordinates 
with respect to the CM. These generators are used to 
build the basis, the interaction, the many-particle kinetic 
energy operator, as well as to evaluate observables. 

We note that the NCSpM, as presented here, is lim- 
ited to interactions that preserve the Sp(3,R) symme- 
try. While this greatly facilitates the use of a group- 
theoretical apparatus and analytical expressions for the 
Hamiltonian matrix elements, which, in turn, makes large 
iV max spaces possible to manage, ab initio calculations 
lie beyond the scope of the current model. The addition 
of symmetry-mixing terms in the interaction is feasible 
and a logical extension of the theory that is under devel- 
opment. Nonetheless, Sp(3, R)-symmetric Hamiltonians 
appear to be particularly suitable to capture the essen- 
tial characteristics of the low-energy nuclear kinematics 
and dynamics. The reason is that such Hamiltonians can 
include the many-particle kinetic energy, X^pf/(2ra), 
the HO potential, rafi 2 r?/2, as well as terms depen- 
dent on the orbital momentum, L = x p^, and 
terms dependent on the mass quadrupole moment, Q2/1 = 
J2i Q2ii(i) = J2i \/ 167r / 5r 2 2 ^2/i(ri), such as the important 
interaction of each particle with the total quadrupole mo- 
ment of the system, \Q ■ Q = \ £\ q(i) • (£\ q(j)) [18]. 

For the purposes of this study, we utilize a novel mi- 
croscopic many-body interaction suitable for large- 7V max 
no-core shell models, which is tied to the long-range ex- 
pansion of the nucleon- nucleon central force V(\ri — rj\) 
(for a derivation, see [22]) and kept as simple as possible 
by considering the most relevant degrees of freedom for 
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a description of highly deformed spatial configurations, 
namely, coupling of single-particle q2 fJj (i) quadrupole HO 
excitations, 

ff 7 = ^^ + ^— J"2 ^ 

given in terms of particle coordinates relative to the CM 1 . 
We take the coupling constant x — Ml/(A^/NfNi) with 
Nfty the total HO quanta of the final (initial) many-body 
basis state. The decrease of x with A/", to a leading order, 
has been shown by Rowe [25] based on self-consistent 
arguments and used in an Sp(3, R)-based study of cluster 
states of 16 O [26] . It is worth noting that such a choice for 
X renders the eigenstates of H 1 ^-independent. Above 
all, the H 1 introduces simple but important many-body 
interactions that enter in a prescribed hierarchical way 
given in powers of a small negative parameter 7, the only 
adjustable parameter in the model. 

The NCSpM, with H 1 , reduces to the established El- 
liott model [19] in the limit of a single valence shell and 
zero 7. In this limit, the model was shown to effectively 
describe rotational features of light nuclei [19 . A suc- 
cessful extension to multiple shells has been achieved and 
applied to the 24 Mg g.st. rotational band [27], where an 
interaction given as a polynomial in Q up through (Q-Q) 2 
was employed. Such an interaction directly ties to our ef- 
fective Hamiltonian 0. In addition, Eq. could be 
understood in terms of a renormalization (e.g., see [28]) of 
the x coupling constant of the \ Y^ij 'Q(j) NN inter- 

action, i.e., f (e^<? - l)/ 7 = i[ X QXo ^f^)]Q^- 
where higher-order terms become quickly negligible for 
a reasonably small 7 and, e.g., for 12 C only terms up 
through k — 1 (g.st.) or k = 3 (Hoyle state) are found 
sufficient. 

Energy spectrum for 12 C and comparison to ab ini- 
tio results. — As the interaction and the model space 
are carefully selected to reflect the most relevant physics, 
the outcome reveals a quite remarkable agreement with 
experiment. The low- lying energy spectrum and eigen- 
states for 12 C were calculated using the NCSpM with H 1 
of Eq. for Ml = 18 MeV given by the empirical es- 
timate « 41/A 1 / 3 = 17.9 MeV. The results are shown 
for iV max = 20, which we found sufficient to yield con- 
vergence. This iV max model space is further reduced by 
selecting the most relevant symplectic irreps (Fig. [I]), 
namely, the spin-zero (S = 0) OMl Op-Oh (04), 2MI 2p- 
2h (6 2), and AMI 4p-4h (12 0) symplectic bandheads to- 
gether with all multiples thereof up through 7V max = 20 



Although a technical detail, it is important to note that Q.Q in 
this expression denotes the Q.Q — (Q.Q) n interaction, where the 
(Q.Q)n term is subtracted from Q.Q to eliminate a spurious shift 
in the zero-point energy by the average contribution, (Q.Q) n , of 
Q.Q within the subspace of n HO excitations |23l I24| . 



of total dimensionality of 4.5 x 10 3 . In comparison to the 
experimental energy spectrum (Fig. [2^t) , the outcome re- 
veals that, for 7 = —1.71 x 10 -4 , the lowest + , 2 + , and 
4 + states of the OhQ Op-Oh (04) symplectic irrep closely 
reproduce the g.st. rotational band; the AMI 4p-4h (12 0) 
irrep successfully describes the Hoyle state and its rota- 
tional structure, while the calculated lowest + of the 
2MI 2p-2h (6 2) is found to lie close to the 10-MeV 0+ 
resonance (third + ) observed in 12 C. We note that our 
model yields a rather compressed spectrum (Figure [2^l, 
dotted-line states) with a small 0^ st -2± spacing as the 
one calculated by early cluster models, which remedy this 
by allowing for alpha-cluster dissociation due to a spin- 
orbit force as discussed in Ref. [30]. However, an overall 
factor of two brings all the NCSpM states under con- 
sideration remarkably close to the experimental values 
(Figure [2^i, solid-line states). This scaling has no effect 
on the underlying physics, as also shown next, because 
an overall factor for H 1 does not change its properties 
and eigenstates, or other associated observables. 

A close similarity is observed when the probability dis- 
tributions for the g.st. rotational band are compared to 
ab initio results when only configurations of zero proton 
and neutron spins (S p , n = 0) are selected (for 0^ st , see 
Fig. [3|i). In particular, NCSpM eigenstates, which are 
^-independent, are compared to symmetry- adapted no- 
core shell model (SA-NCSM) calculations [13] with bare 
JISP16 realistic interaction for Ml = 20 MeV (around 
the minimum of the calculated binding energy for 12 C) 
and an iV max = 8 model space. This space appears to 
be sufficient to yield convergent results for the g.st. rota- 
tional band for both models. The close agreement shows 
that among all possible S PjIl = configurations present 
in the SA-NCSM, only the states of the (0 4) slice used 
here appear dominant. It also points to the fact that the 
schematic interaction used in NCSpM has effectively cap- 
tured most of the underlying physics of the realistic in- 
teraction important to the low-energy nuclear dynamics. 
Indeed, for the same iV max and Ml, NCSpM observables, 
such as g.st. matter rms radius, B(E2; 2+ 0+ s J, and 
Q 2 +, reproduce the ab initio counterparts as much as 
91%, 83%, and 80%, respectively. 

For the Hoyle state and its rotational band, larger 
spaces are needed, e.g., nonnegligible configurations ex- 
tend to iV max = 18 (Fig. |3Jd), which is within a reach 
of next-generation NCSM models. While the predomi- 
nant component of the lowest + state is at 0h£l (n = 0) 
and manifests an evident oblate shape, the Hoyle state 
peaks around SMI (n = 8) with a clear evidence for a 
prolate shape deformation with (16 0) being the largest 
contribution. Moreover, this state emerges from a 4p-4h 
shell-model configuration, that is, the (12 0) bandhead 
(Fig. [l] set I) is realized by an alpha-particle configura- 
tion - spatially spherical (0 0), spin zero, and isospin zero 
- in each of the three lowest HO shells (implying spatial 
displacement). This together with the strong prolate de- 
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FIG. 3: Probability distribution for 12 C as a function of the 
n total excitations, (a) of the lowest + state as calculated 
by the NCSpM (left) and the SA-NCSM (right), and (b) of 
the Hoyle state as calculated by the NCSpM. The dominant 
shape deformations, specified by (A/i), are also shown. 
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FIG. 4: Dependence of the 12 C NCSpM energy spectrum on 
(a) the model space (A^ max ) for 7 = — 1.71 X 10 -4 and (b) the 
7 model parameter for iY max = 20. 



formation of this state (see also Q of its 2 + excitation in 
the next section) supports an underlying alpha-particle 
cluster structure of this + state and its rotational band. 

It is interesting to point out that the AhVt 4p-4h (12 0) 
places the controversial 2 + excitation within the lower 
region of the 9.6-12 MeV range of observations [2 - 6 . For 
comparison, recent ab initio iV max = 8 NCSM calcula- 
tions, while achieving a remarkable reproduction of the 
g.st. rotational band, yield the second + and 2 + states 
around 13 MeV and 15 MeV, respectively [10 , thus be- 
lieved not to be associated with the Hoyle state but with 
higher-lying states of that spin-parity. Indeed, consistent 
with current ab initio observations, the NCSpM outcome 
demonstrates a large sensitivity of the energy of the Hoyle 
state and its 2 + excitation on the model space (Fig. Bp). 

Furthermore, a study of the effect of 7 on the ^C 
energy spectrum (Fig. ^p) reveals that the additional 
degree of freedom associated with the 7 model parameter 
is in fact substantially limited by the lowest + states 
(with only a small effect on the g.st. rotational band). 
Indeed, given the dramatic variation with 7 for the 0^ 
and O3" (Fig. [4)3, insert), there is only a small window of 
reasonable 7 values where observables are also found in 
agreement with experiment, as shown in the next section. 

Finally, to show the important role of the 2h£t and 
Ahfl irreps considered in this study, we include additional 
5 = symplectic irreps into the model space. In doing 
this, we indeed find no additional low- lying states for the 
7 parameter chosen here. For example, for 2h£t (2 4), Ah£t 
(8 2), (44), and (0 6), 6HQ (14 0), as well as 8ftfi (16 0), 
the lowest state (0 + ) lies much higher than 18 MeV. 



Additional observables. — The model with the se- 
lected 7 = —1.71 x 10 -4 successfully reproduces other 
observables for 12 C that are informative of the state 
structure, such as r rms point-particle matter rms radii, 
Q electric quadrupole moments (Table [i]) and B(E2) 
transition strengths (Fig. [2^t) [in general, for 7 from 
—0.01 to — 10 -5 , e.g., r rms for 0^ st and 0^ increases 1.4 
times, as well as B(E2; 2f — > 0~^ st ) and Q 2 + increase 
four times]. However, nonzero B{E2; Oj 2+) and 
M(E0; O2 —>0t) strengths can only result from mixing 
of Sp(3,R) irreps, which requires an Sp(3,IR) symmetry- 
breaking interaction. But this can enter perturbatively. 
Indeed, a useful information on these quantities can be 
derived from a toy model that mixes the (12 0) irrep 
into the (04) irrep of the g.st. and that reveals that less 
than a 2% mixing can already yield the correct order of 
magnitude, namely, B(E2; 0^ -» 2^) = 8.7 W.u. and 
M(E0; O2 0^) = 2.04 efm 2 as compared to experi- 
ment, 8.0(11) W.u. and 5.4(2) efm 2 , respectively. 

We also find matter rms radius for the ground state 
(Table [i]) close to the value deduced from experiment 
in Ref. [3T] , Interestingly, our calculations yield mat- 
ter r rms = 2.93 fm for the Hoyle state. While this re- 
sult drastically differs from predictions of cluster models 
[3 [34j [35], it is close to a recent value deduced from 
the experiment, 2.89(4) fm [33]. The NCSpM outcome 
yields a 2^ electric quadrupole moment very close to the 
experimental value, and a large negative one for the 2 + 
above the Hoyle state (Table [i]) indicating a substantial 
prolate deformation for the Hoyle and 2 + states. Such a 



TABLE I: Point-particle rms matter radii (fm) and electric 
quadrupole moments (efm 2 ) for 12 C and 24 Mg. Experimental 
data taken from: a Ref. [31., b Ref. [32], c Ref. [33], d Ref. [29]. 
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prolate deformation, albeit not so pronounced, has been 
also suggested by the ab initio lattice EFT [TT] . 
NCSpM applications to other nuclei. — The out- 
come of the present analysis is not limited to 12 C. The 
model we find is also applicable to low-lying states of 
other light nuclei, and, in particular, we show here the 
g.st. rotational band and low-lying + states of 24 Mg 
(Fig. J2Jd). These states together with associated observ- 
ables (see also Table [i]) are well described by the NCSpM 
with 7 = —0.67 x 10 -4 in an 7V max = 24 model space us- 
ing only a few spin-zero symplectic irreps, namely, OhQ 
Op-Oh (8 4) for the g.st. rotational band, OMl (4 6), (6 2), 
and (0 8), 2KI 2p-2h (14 2), as well as 4fifi 4p-4h (20 0). 
This points to the efficacy and extensibility of the fully 
microscopic NCSpM model, which has indeed captured 
most of the important physics. 

In short, we carried forward a no-core shell- model 
study with a schematic many-nucleon interaction to fur- 
ther unveil the underlying physics behind various phe- 
nomena important to the low-energy nuclear dynamics 
of 12 C. We showed, for the first time, how both collec- 
tive states and states suggested to have cluster-like sub- 
structures emerge out of a fully microscopic, shell-model 
framework, thereby providing a novel and essential per- 
spective on the controversial Hoyle state. 
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